Stem cell population asymmetry can reduce rate of replicative aging 
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Cycling tissues such as the intestinal epithelium, germ line, and hair follicles, require a constant 
flux of differentiated cells. These tissues are maintained by a population of stem cells, which generate 
differentiated progenies and self-renew. Asymmetric division of each stem cell into one stem cell and 
one differentiated cell can accomplish both tasks. However, in mammalian cycling tissues, some stem 
cells divide symmetrically into two differentiated cells and are replaced by a neighbor that divides 
symmetrically into two stem cells. Besides this heterogeneity in fate (population asymmetry), stem 
cells also exhibit heterogenous proliferation-rates; in the long run, however, all stem cells proliferate 
at the same average rate (equipotency). We construct and simulate a mathematical model based on 
these experimental observations. We show that the complex steady-state dynamics of population- 
asymmetric stem cells reduces the rate of replicative aging of the tissue -potentially lowering the 
incidence of somatic mutations and genetics diseases such as cancer. Essentially, slow-dividing stem 
cells proliferate and purge the population of the fast-dividing -older- cells which had undertaken 
the majority of the tissue-generation burden. As the number of slow-dividing cells grows, their 
proliferation-rate increases, eventually turning them into fast-dividers, which are themselves replaced 
by newly emerging slow-dividers. Going beyond current experiments, we propose a mechanism for 
equipotency that can potentially halve the rate of replicative aging. Our results highlight the 
importance of a population-level understanding of stem cells, and may explain the prevalence of 
population asymmetry in a wide variety of cycling tissues. 



I. INTRODUCTION 

As cells replicate they accumulate damage, such as so- 
matic mutations, which can eventually manifest as a mal- 
ady; cancer is the deadliest culmination of such replica- 
tive misfortunes [iHS]. Cells that undergo many rounds 
of divisions are more likely to develop cancer; particularly 
susceptible are the somatic stem cells of tissues that con- 
stantly cycle, i.e. the mammalian intestinal epithelium, 
epidermis, and germ line [6l[7]. 

Stem cells are defined for their capability to generate 
more stem cells -self-renewal- and daughter cells that 
differentiate [8] . Asymmetric divisions of a stem cell into 
one stem cell and one differentiated cell will satisfy both 
these objectives: replenish the stem cell pool and gener- 
ate the differentiated progenies cycling through the tis- 
sue [9l llOj . However, lineage-tracing experiments have 
demonstrated that in a diverse set of mammalian cy- 
cling tissues, stem cells divide symmetrically, generating 
daughter cells that acquire the same fate [TTMlSj . The 
balance between proliferation and differentiation is main- 
tained at the population level: some stem cells divide 
into two stem cells, whereas others generate two differ- 
entiated cells and are lost [H]. To maintain homeostasis 
-constancy in cell number and tissue organization, both 
outcomes occur with frequency one-half. 

The justification for the prevalence of population 
asymmetry in cycling tissues remains elusive. The as- 
sertion that symmetric divisions are required for regen- 
erating the tissue following injury and during develop- 
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ment [101 E! i motivates the capability but not the con- 
tinual use of this strategy in homeostasis. For instance, 
although germ line stem cells in the Drosophila ovaries 
are capable of symmetric division when a stem cell is lost, 
during normal proliferation all divisions are asymmetric 
[TPj . We propose that population asymmetry reduces the 
rate of replicative aging in mammalian cycling tissue. 

Herein, we formulate a model and simulate the com- 
plex dynamics of stem cells with population asymmetry 
and external regulation: stem cells lost to differentiation 
are replaced by their neighbors. The source of stochas- 
ticity in stem cell fate is the asynchronous waiting time 
before division. We show that heterogenous stem cells - 
cells that cycle at different rates- with population asym- 
metry undergo fewer divisions for the same flux of differ- 
entiated cells compared with stem cells that use division 
asymmetry. 

The degree of heterogeneity does not need to be large: 
the slowest dividers are only four times slower than the 
fastest dividers in the population, consistent with exper- 
imental observations. A single stem cell's cycling-rate 
fluctuates, but its long-time average is equal to the pop- 
ulation average. All stem cells are therefore equipotent; 
no single cell has a proliferative advantage over the others 
in the long run. 

Going beyond current experiments, we show that the 
replicative-aging rate can be further reduced using the 
following strategy: stem cells by default divide faster 
than the average rate of cycling in the population. Ex- 
ternal stimuli slow down the division rate, a trait which 
is encoded in methylation patterns and passed on to the 
daughter cells. After a random number of divisions the 
methylation signal is lost and the daughter stem cells re- 
vert to their default behavior. This simple strategy can 
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FIG. 1. Heterogeneous proliferation rates. A) Two stem cells 
(orange) divide asymmetrically three times at a constant rate 
generating six differentiated progenies (blue). Both stem cells 
will have aged by the same amount (three divisions). The 
phylogeny tree of each stem cell -neglecting the branches cor- 
responding to somatic cells- is simply a straight line. B) The 
stem cell on the left undergoes four asymmetric divisions fol- 
lowed by a symmetric division into two differentiated cells. 
Stem cell on the right divides symmetrically into two stem 
cells at a lower rate. The net product is the same as A. The 
two stem cells, however, have aged by only one division. With 
symmetric divisions, the phylogeny of a stem cell is more com- 
plicated. The highlighted part of the tree depicts a stem cell 
that has undergone 8 divisions and is younger than some of 
the other descendants. 




Top-view Side-view 

FIG. 2. Population asymmetry. Stem cells form an epithe- 
lium attached to the basal membrane (side-view and top- view 
schematic shown here). A stem cell divides symmetrically 
into two differentiated cells and leave the basal layer. The 
lost stem cell is replaced by one of its neighbors which divides 
laterally into two stem cells. Both types of symmetric division 
occur with frequency 1/2. The net product of the two divi- 
sions is two differentiated cells and two stem cells. Above, 
the third stem cell from the left has divided symmetrically 
into two differentiated cells and is leaving the basal layer. Its 
neighbor to the right is dividing symmetrically into two stem 
cells to replace the lost stem cell. 



reduce the rate of replicative aging by a factor of two; we 
conjecture that this mechanism is used in mammalian 
cycling tissues. 

II. RESULTS 

We will refer to replicative aging -the cumulative num- 
ber of divisions undergone by a cell- simply as 'aging'. 
Cell aging without replication can occur through post- 
mitotic mutations in mitrochondrial DNA f20], stress- 
induced DNA damage, or accumulation of altered pro- 
teins pT|[^. 

If a population of stem cells divides synchronously, all 
the stem cells will age by the same amount regardless of 
their ancestry; both division asymmetric and population 
asymmetric strategies are equivalent. A slight degree of 
heterogeneity in the division rates, however, can lower the 
average aging rate, if symmetric divisions are allowed, 
where some stems cells can replace others (Fig.l). In 
such a population, the fast-dividing cells take over most 
of the proliferative burden and are eventually replaced 
by slow-dividing cells. 

We consider an idealized model of cycling tissues 
(Methods) : a fixed number of stem cells form an epithe- 
lial basal layer [23) . After a random waiting-time, stem 
cells divide asynchronously -along their out of plane axis- 
into two differentiated progenies, which leave the layer. 
In response, a randomly-chosen neighbor symmetrically 
divides in the lateral direction and replaces the lost stem 



cell -externally regulated division. This model closely 
resembles the actual dynamics of mammalian cycling tis- 
sues [T21 - [T5] . Since a division into two differentiated cells 
is always accompanied by another into two stem cells, 
the frequency of both types of division is 0.5, ensuring 
a constant flux of differentiated cells and self-renewal of 
the stem cell population. 

The above proliferation strategy results in dynamics 
that closely resemble the 'stepping-stone' model of neu- 
tral drift in population dynamics [121 [21] (see also 'voter' 
models in statistical physics [25]). If we initially label all 
cells in the population with a different color (Fig. 3), over- 
time, we observe that most cells leave the stem cell layer, 
and are replaced by their neighbor. Those that remain, 
grow in number, forming clusters of daughter cells (of the 
same color) called clones. If all stem cells are equivalent 
and cycling at the same rate, no clone is more likely to 
survive than its neighbors. The boundaries between adja- 
cent clones diffuse back and forth, roughly expanding as 
\/t\ the size of the surviving clones therefore increases as 
t (for a rigorous derivation and the log correction to this 
scaling, see EQ). Eventually, in a finite population, one 
clone will take over, achieving fixation; in Fig. 3, a stem 
cell initially labeled red fixates. The resulting distribu- 
tion of clone sizes is scale-free: the distribution of the 
ratio of clone sizes to the average clone size at any given 
time is constant. Statistical signatures of clone sizes have 
been used to experimentally confirm the above dynamics 
in cycling tissues in vivo [16J. 
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FIG. 3. Clonal expansion and loss of heterogeneity. Cells 
are initially colored at random and their evolution followed. 
Eventually the clone of a red stem cell fixates. The time steps 
from left to right and bottom to top are, t=l, 20, 60, 100,140, 
200, and 260. The stem cells have heterogenous division rates, 
a = 0.15. 



Heterogeneity 

Next, consider a heterogenous population, where stem 
cells cycle at different rates: stem cell i divides during 
each time step 5t = 1 at rate drawn from a normal 
distribution with mean 0.5 and standard deviation a. 
The mean cycling-rate is held fixed to ensure a constant 
flux of one differentiated cell per stem cell per unit of 
time 5t (Methods). The clones are no longer equivalent; 
some are descendants of fast-dividing cells and others of 
slow-dividing cells. Clusters of fast-dividing cells con- 
tribute more to the flux since they cycle more quickly. 
However, at their interfaces, the slow-dividing clones in- 
vade the fast-dividing clones, since at any time-step, the 
fast-dividing cells are more likely to undergo terminal 
differentiation and be replaced by a neighbor. The fast- 
dividing cells that are contributing more to the flux -and 
thereby aging faster- are systematically removed from 
the population and replaced by younger cells. 

Define the population-averaged aging rate p{t) as the 
ratio of the total number of divisions starting from time 
t = undergone by the stem cells present in the pop- 
ulation at time t to the total number of somatic cells 
produced since t = Q. p — 1 for homogenous population 
-where all cells proliferate at the same rate, since all the 
stem cells contribute to the flux equally and on average 
one differentiated cell is produced by each division. Sim- 
ilarly, p = 1 for a population with division asymmetry. 
Surprisingly, our simulations show that the aging rate at 
early times is lower than that of a homogenous or division 
asymmetric population, p < 1 (Fig.4A). 

We can understand the eventual reversion to p = 1 
by considering the dynamics more carefully. The ini- 
tially heterogenous population of cells flrst gives rise to 
an abundance of slow-dividing cell, since they are more 
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FIG. 4. Heterogeneous but not equipotent. A) Population- 
averaged aging rate p as a function of time. Initially p < 1 
before reverting to its asymptotic value of 1. 50 runs are su- 
perimposed for each value of a. The thick colored lines are 
the average over the runs. Solid-blue lines are the analytical 
results of the phenomenological model. B) Frequency of ob- 
serving a cell that will divide with probability r in time St for 
the first 120 time steps superimposed {a = 0.15). The distri- 
bution narrows over time, eventually becoming a (5-function 
as one clone fixates. C) Clone sizes in the population as a 
function of time for the same parameter values. Some clones 
expands as others shrink. D) Average clone size as a function 
of time. There is a significant deviation at early times from 
the linear prediction of the neutral drift model. 



likely to invade their neighbors. As the number of slow- 
dividing cells grow -clonal expansion, their cycling-rate 
increases to ensure a constant flux of differentiated proge- 
nies. The initially Gaussian distribution of cycling-rates 
in the population develops a spike at the slowly-dividing 
rates, the spike grows and reverts to the mean (Fig.4B). 
Eventually the distribution becomes a (5-function at 0.5 
when flxation has occurred and all the stem cells in the 
population have the same proliferation rate. 

Phenomenological model 

A simple phenomenological model can elucidate the es- 
sential mechanisms for lowering the rate of aging. We 
include the following key facts in the model: 1) The 
boundaries of the slow-dividing clones grow faster than 
would be expected from neutral dynamics (see Fig.4D). 
At early times, the boundaries grow on average with ve- 
locity 0.5 — r, where r is the cycling-rate. 2) As the num- 
ber of slow-dividing cells increases, their cycling-rate also 
increases, approaching the population-average rate, 0.5. 
We assume that fractional rate of increase of the cycling- 
rate r is the same as the fractional increase in the number 
of cells dividing at rate r. 
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Using a Lagrangian framework, we follow the evolution 
of the cells that at i = cycled at rate tq; at time t, 
the number of these cells in the population is denoted 
as N{ro,t) and their cycling-rate as r{ro,t). In 2D the 
dynamical evolution is given by, 



1 



dN{ro,t) 



2^Niro,t) dt 



(1) 



dr{ro,t) 1 9iV(ro,0 



U.t) (2) 



dt N{ro,t) at 

The number of stem cells cycling at rate r' at time t is 

AA(r',i)= [ N{ro,t)-^S{r'-riro,t))dr,. (3) 

Eqs. [T]and[2] can be combined into the following non- 
linear ODE -denoting N{rQ,t) as N for simplicity. 



cPN /I 1 \dN 1 fdN 



dt"^ \2N Jn I dt N\dt 



(4) 



The initial distribution of proliferation rates in the sim- 
ulations is a Gaussian with mean 0.5 and standard de- 
viation a. For simplicity, we follow two population of 
cells from this distribution: fast-dividing cells cycling at 
rate r/ — 0.5 + cr, and slow-dividing cells cycling at rate 

= 0.5 — cr. Simulations (Fig. 4) show a monotonous de- 
crease in the second moment of the cycling-rate distribu- 
tion -the mean is fixed by construct; the third moment, 
however, peaks when p is at its minimum. We claim 
that to the lowest order, the rate of aging at time i is a 
function of the third moment of the distribution of the 
cycling-rates in the population at time t. The more lop- 
sided the distribution, the more prevalent the presence 
of slow-dividing clusters that invade their fast-dividing 
neighbors, ensuring p < 1. The analytical estimate of 
p{t) is 



l-p{t)^[N{T,,t){]^~r{rs,t)f-N{rf,t){T{Tf,t)~\fY 

(5) 

The exact shape of the curve p{t) depends on system 
size. The loss of heterogeneity -rate of reversion of p{t) 
to its asymptotic value- will be slower for larger systems 
since there is a larger pool of initially slow-dividing cells. 
The analytical theory is a good fit to the simulations 
results (Fig.4A , Fig.7C); only two fitting parameters are 
used (see Methods). 



Heterogeneous and equipotent 

To sustain the low rate of aging, we need to ensure that 
heterogeneity is not lost by clonal expansion, and some- 
how reintroduced into the system. Inspired by methy- 
lation as a method of passing epigenetic information to 



daughter cells (see Discussion), we assume that at each 
division there is a probability that the inherited methy- 
lation patterns are lost and the cell's cycling-rate ran- 
domized (see Methods for implementation). The average 
number of divisions (time-scale) to losing inherited infor- 
mation is denoted as r^q. If we average the behavior of 
the cells over time-scales much longer than Tg^, all cells 
are equipotent: there are no slow or fast dividers, all 
cells proliferate at the average rate. With equipotency, p 
reaches a steady-state value poo, which is approximately 
equal to p(Tgg) in the heterogenous but non-cquipotcnt 
case (Fig.5A). 

A steady-state rate of aging poo < 1 does not imply 
that there is a steady-state distribution of cycling-rates, 
i.e. a constant number of slow and fast dividers. The 
dynamics is more complicated: as in the non-equipotent 
case, slow dividers expand as their cycling-rate increases 
(Fig.5B). Equipotency, however, ensures that somewhere 
else in the population, a new slow-divider emerges and 
expands; the cycle continues (Fig. 6). Time-averaging the 
cycling-rate distribution for a sufficiently-long periods, 
yields a static skewed distribution -non-zero third mo- 
ment (Fig.5C). There is an abundance of slow-dividers 
(tail of the distribution) which are compensated by an 
abundance of cells that divide only slightly faster than 
the average. 

Equipotency also restores the scale-free dynamics; 
stem cells behave like a homogenous population for time 
scales larger than Tgg. For example, the scale-free charac- 
ter of the clone size distributions, and the linear growth 
of the average clone size is restored (Fig.5D), consistent 
with experimental observations [ TTHTB] . 



Optimizing 

If proliferation-rates randomize too quickly, slow- 
dividing cells do not have sufficient time to replace their 
faster dividing neighbors. If the randomization rate is too 
slow, the slow-dividing clones expand too much, reduc- 
ing the heterogeneity, and increasing the rate of aging. 
There is an intermediate optimal time-scale for equipo- 
tency. We examined the dependence of the steady-state 
value of aging rate poo on Teq (Fig.7A). The optimal value 
of T*g that minimizes aging rate is approximately the 
time required to achieve the transient minimum p in a 
heterogenous non-equipotent population; t*^ is the point 

of extremum in Eqjij ^^^^|r*^~o- 

To further optimize p, we consider other feasible mech- 
anisms: 

Bimodal distribution of cycling-rates. Thus far, 
equipotency was implemented by randomizing cycling- 
rates after a few divisions. Can a different mechanism do 
better? In Fig.7B, cells establish equipotency by reset- 
ting their cycling- rate to above-average tr -I- 0.5. Physi- 
cally, this corresponds to stem cells that are programmed 
by default as fast-dividers. However, global signaling es- 
tablishes a methylation pattern where most cells divide 
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Cycling Rate Time 

FIG. 5. Heterogenous and equipotent. A) Aging rate p as 
a function of time for a = 0.15 and four different values of 
Teq {Ne ~ 1, see methods). 50 runs shown and averaged. B) 
Distribution of cycling-rates superimposed for 20 time steps 
(3980 to 4000) for r^q = 250. At steady states waves of slow 
dividers are generated and sent towards the mean, where they 
annihilate (see Fig. 6). C) Distribution of proliferation-rates 
averaged for 1000 steps (3000 to 4000). The time-averaged 
distribution looks Gaussian but is asymmetric (non-zero third 
moment). The inset shows the ratio of the actual distribution 
to a Gaussian fit. There is an excess of slow divider cells -at 
the tail on the left-side. D) Average clone size as a function 
of time for Teq — 1. The scale- free behavior of the neutral 
drift models is recovered when division rates are randomized. 




FIG. 6. Maintaining heterogeneity. The color-coding cor- 
responds to the division probability of each cell in time St 
(cycling-rate). Each cell randomizes its division rate (cr = 
0.15) with Teq = 50. A slow-divider in the middle proliferates 
while slowing down, eventually reaching the mean prolifera- 
tion rate 0.5. A new cluster of slow-dividing cells has formed 
in the bottom-left corner and is expanding. The time steps 
shown are from left to right and bottom to top, t=90 to 140 
in increments of 10. 



slower; after a random number of divisions the cell loses 
the methylation program and reverts to its default be- 
havior. 

This approach results in a bimodal distribution of 
cycling-rates: slow-dividing cells that replace their neigh- 
bors and proliferate, balanced by the fast-dividing cells 
that are constantly injected into the population when 
cells reset their cycling-rates. For optimal equipotency 
time-scale, r*^ and a factor of 6 difference between the 
proliferation-rates of the fastest and slowest dividers, we 
find Poo =0.4. This, implies approximately a two-order 
of magnitude (99%) reduction in incidence of five so- 
matic mutations -on average required for onset of cancer- 
compared with homogenous or division-asymmetric pop- 
ulations (Methods). For our choice of parameters, the 
optimal time-scale for equipotency is now much shorter 
(around 3-5 divisions). Quick addition of fast-dividers is 
required to balance the growing number of slow-dividing 
stem cells. This strategy is so effective that we conjec- 
ture that it is potentially implemented in cycling tissues 
(see Discussion). 

Biassed replacement. We consider the possibility 
that beside terminal differentiation a stem cell can also 
proliferate by forcing the terminal differentiation of a 
neighbor. Fig.7C shows the aging-rate where a stem cell 
cycling at rate r, proliferates with probability 1 — r. A 
modest improvement is observed (4% for our choice of 
parameters) since the fast-dividing cells are even more 
likely to be purged from the population. 

One-dimension. Finally, we consider a population of 
stem cells constrained to one-dimension. Although the 
steady-state aging rate poo is not significantly different 
from that of the two-dimensional case, the dynamics of 
clonal expansion and loss of heterogeneity is qualitatively 
different (Fig.7D). For little or no heterogeneity, the clone 
boundaries expand by diffusion: the average surviving 
clone size increases as ^/t (as opposed to t in 2D). As 
heterogeneity increases the boundaries drift rather than 
diffuse and the clone size scaling becomes closer to linear. 
The slower clonal expansion implies that heterogeneity is 
lost at lower a rate, which maximizes the transient period 
where p < 1. ID systems can help optimize the aging rate 
for finite durations for systems without equipotency. 



III. DISCUSSION 

Wc found that population asymmetry in stem cells 
with heterogeneous division rates and external regula- 
tion can substantially reduce the cumulative replicative 
aging of the tissue; potentially lowering incidence of so- 
matic mutations and proliferative-maladies such as can- 
cer. Our analysis revealed a non-trivial steady-state dy- 
namics: slow-dividing cells constantly replace their fast- 
dividing (older) neighbors who had shouldered the brunt 
of the proliferative demand. A key requirement was 
equipotent stem cells: in the long run no cell had a pro- 
liferative advantage over the others; every cell divided 
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FIG. 7. Minimizing rate of aging. A) Asymptotic value of aging-rate p as a function of T^q for different degrees of heterogeneity 
a {Ne — 10; Methods). Optimal r^q is at some intermediate value. B) Instead of randomizing, the cell's cycling-rate is reset 
to a + 0.5. (Top) p as a function of t for different values of r^g; a = 0.15. Minimum p occurs at the high equipotency rate of 
Teq = 2.5. The bimodal distribution of the division rates is shown on the bottom -every 5th step for the first 2000 iterations: 
cells are either fast dividers or slow dividers. The inset shows the spatial heterogeneity of the two types of stem cells (color- 
coding as in Fig. 6). On the bottom-left, we plot poo as a function of r^q for various value of a. C) Decision of fate is biassed so 
that slow-dividing cells are more likely to replace a neighbor. (Top) p as a function of time and no equipotency. The solid blue 
lines are the fits from the phenomenological model (Methods). (Bottom) Equipotency for four different values of Teq; cr — 0.15. 
D) ID system of size 10000x1. (Top) p as a function of time. Loss of heterogeneity (clonal expansion) is slower in ID. (Bottom) 
Average clone size as a function of time for two values of a. The clones grow sub-linearly. B) Top, C), and D) Top, multiple 
trajectories are 50 runs and their average (thick line). 



at the same average rate. Equipotency was enforced by 
assuming that after a random number of divisions, the 
resulting daughter cells lost their inherited traits. Go- 
ing beyond current experimental findings, we conjectured 
that if equipotency resets the cycling-rate of the cells to a 
default value that is larger than the mean division-rate of 
the population then a significant reduction in the aging 
rate can be achieved. 

Stem cell population asymmetry using external regu- 
lation -where terminal differentiation of a stem cell trig- 
gers the expansion of a neighbor- has been verified using 
lineage-tracing experiments in the mammalian intestinal 
crypt and germ line |12H15j . These quantitative studies 
have demonstrated that these tissues are maintained by 
an equipotent quickly-cycling population of stem cells. 
Do these populations also exhibit heterogeneity? Direct 
observation of heterogeneity in the intestinal epithelium 
was made by Pruitt et al. |27j who using tamoxifen in- 
duced Cre-recombinase activity on the Mcm2 gene ob- 
served cell-cycle times ranging from 1 to 4 days, consis- 
tent with the degree of heterogeneity assumed above. 

Despite a lack of direct evidence of heterogeneity 
of cycling-times in the germ line, there are observa- 
tions of heterogeneous gene-expression in undifferenti- 
ated cells of the mammalian spermatogonia |28j . Two 
types of undifferentiated cells and Apr show an ex- 



pression level low in NgnS and high in GFRal, whereas 
undifferentiated Aai cells have the opposite expression 
levels. Aai cells are thought to be primed for dif- 
ferentiation and could correspond to the fast-cycling 
cells in our model. Short-term lineage tracing experi- 
ments have shown that the expression patterns can re- 
versibly change with GFRaV-° / N gnZ^^ cells becoming 
GFRal'''' /Ngn3'°, ensuring equipotency [H]. When As 
and Apr cells commit to differentiation, they are replaced 
by fragmentation of a neighboring Aai cell that reverses 
its expression patterns p4] . 

Intestinal crypt stem cells are also believed to exhibit 
dynamics where the terminal differentiation and loss of 
one stem cell results in expansion of a neighbor. Wnt sig- 
naling 1221 [3D] -in particular a global Wnt gradient along 
the crypt [3T]- is thought to control the average prolifer- 
ation rae[32]; as do insulin-like growth factors like IGFl 
|27) . Nearest-neighbor communications seem to be me- 
diated by Notch signaling and lateral inhibition. In the 
crypt, stem cells laterally inhibit their neighbors to en- 
sure that the neighbor's fate at division is differentiation 
|33) . A biological mechanisms for the biassed expansion 
proposed above (Fig.7C) might be the build-up of excess 
Notch ligands by slow-dividing cells. 

We assumed that cycling-rates were passed on to 
daughter cells after division. Moreover, to establish 
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equipotency, we required that the rates randomize af- 
ter a few divisions; the inherited information forgotten. 
Both of these requirements can be implemented through 
epigenetic methylation tags [331I3S]. Yatabe et ah [3B] 
studied methylation patterns in the human colon crypts 
and observed somatic inheritance, random methylation 
changes, and evidence of population asymmetry. 

Our model is an idealized picture of stem cell prolif- 
eration. We only focused on aging in the somatic stem 
cells and not the differentiated progenies as they are con- 
tinually discarded. We also assumed that stem cells di- 
vided symmetrically into two differentiated cells. The 
product of this division can also be thought of as transit- 
amplifying progenies that undergo few more rounds of 
division before terminal differentiation. Inclusion of this 
detail will not change our results. 

In our model, the only source of stochasticity of stem 
cell fate is the random waiting-time prior to division. On 
division, the fate is always terminal differentiation and 
loss; a stem cell proliferates only if its neighbor is lost. A 
more realistic model will introduce a degree of stochas- 
ticity to the fate itself (as was done for the biassed in- 
vaders) . A stochastic fate with external regulation seems 
to be more complex than the deterministic fate consid- 
ered here. Moreover, aging-rate will still be reduced if 
the probability of terminal differentiation is larger than 
that of proliferation. 

The geometry of a 2D and ID epithelium of stem cells 
used in the model is an over-simplification of real tissues. 
Although stem cells seems to be predominantly confined 
to a thin -few cells deep- layer in cycling tissues, they 
are usually interspersed with other types of somatic cells. 
For example, stem cells at the base of the crypt (Lgr5 
cells) are interspersed by Panth cells, and in the germ 
line by large somatic Sertoli cells '15]. A more compre- 
hensive model should account for this added complexity, 
however, the results will not be qualitatively different. 
Previous models of population asymmetry that have re- 
lied on effective simplified 'lattices' of stem cells are in 
good agreement with experiments |16j . 

Notion of heterogeneity can be employed in a differ- 
ent context. Most populations of stem cells contain qui- 
escent cells [371 131] that cycle so slowly that are effec- 
tively dormant. Despite suggestions that quiescent stem 
cells can increase longevity by replacing actively-dividing 
stem cells [35] , lineage-tracing experiments [THl US] sug- 
gest that these cells play little to no role in homeostatic 
cycling of tissues. Their role seems to be more prominent 
in tissue repair following a catastrophic perturbation [18] 
and during development j37j. The heterogeneity used in 
our model is subtle and the slow-dividing cells are in no 
way quiescent. To ensure that extremely slow-dividing 
cells were not essential for our results, we repeated the 
simulations while constraining the cycling-rates to ±3ct 
and observed the same reduction in the aging rate. For 
instance, constraining the division probability of the bi- 
modal strategy (Fig.7B) to 0.2-0.8, resulted in a popula- 
tion where the fastest dividers proliferated only 3 times 



faster than the slowest dividers; the dynamics still gen- 
erated p = 0.87, a potential 50% reduction in incidence 
of cancer. 

Our results motivate a population level understanding 
of stem cells. Because of heterogeneity and the complex 
dynamics, static gene expression profiles provide only 
limited insight into stem cell phenotypes. Study of popu- 
lation epigenetics and population dynamics of stem cells 
can bring fresh insight into the delicate interplay that 
occurs in homeostatic cycling tissues. 



IV. METHODS 
Simulations 

An asynchronous cellular automata algorithm was 
used: A 2D hexagonal lattice of N points with periodic 
boundaries is generated and perturbed by changing the 
position of each point in x and y directions by addition 
of a random number drawn from a Gaussian distribution 
with zero mean and standard deviation 30% of the lat- 
tice constant. The A'' cells are constructed by a Voronoi 
decomposition of the randomized lattice; two cells that 
share a common edge are denoted as neighbors. Each 
cell is initially assigned a division rate r -the probabil- 
ity that the cell will divide in time-step 5t = 1- drawn 
from a random distribution with mean 0.5 and standard 
deviation a. 

In each time step, the following sequence is executed 
N times: 1) A cell is chosen at random from the lat- 
tice. With probability r (the cell's division probability) 
it divides terminally into two differentiated cells. 2) If 
the division occurs, one of its neighbors is chosen at ran- 
dom to divide into two stem cells and replace the lost 
cell. The daughter cells inherit the division rate, and 
the age (number of prior divisions) of their parent plus 
one. 3) If equipotency is turned on, at every division 
with probability Nf./Teq -where Teq is the time-scale for 
equipotency and N^, the number of steps required prior 
to equipotency- the cell advances one step to resetting 
its division rate. After completing N,, steps, the division 
rate of the cell is drawn from random Gaussian distribu- 
tion with mean 0.5 and standard deviation a. The new 
division rate is inherited by all daughter cells in the sub- 
sequent divisions. At the end of each time step, every 
cell's division rate is modified by the same amount to en- 
sure that the population-average of r is 0.5. On average, 
at every time step each a total of N differentiated cells 
are generated. 

Above algorithm was modified for two cases: 1) In 
Fig.7B (Optimization) equipotency is introduced by al- 
ways resetting the cell's division rate to 0.5 -I- a. More- 
over, the division rates are constrained to the range ±3(t; 
2) In Fig.7C (Optimization) the selected cell divides sym- 
metrically into differentiated cells (step 1 above) with 
probability r; with probability 1 — r it divides into two 
stem cells and removes a random neighbor. 3) In Fig.7D 
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(Optimization) the above algorithm was modified to run 
for a one-dimensional lattice. Nf, = 1 for all plots with 
equipotency, except for Fig.7A, where N^. = 10. The al- 
gorithm was implemented in Matlab R2011b. We used a 
100x100 and 10000x1 lattice for the 2D and ID simula- 
tions respectively (N = 10000). For generating the mo- 
saic plots of the lattices (Fig. 3, 5, inset in 7B) a 30x30 
lattice (N=900) was used. 



Estimating incidence of cancer 

Cells need to accumulate a series of independent muta- 
tions for tumorigenesis ^3j. Assume each step (mutation) 
i is an independent Poisson process with rate ki. The 
probability that step i will occur in the time interval t 
to t + dt is an exponential distribution, ^^(i) = kie~'^^^ 
-known as dwell time. The probability that step 2 oc- 
curs at time t is contingent on step 1 having occurred at 
some time before t' < t, P2{t) = /q Wi{t')w2{t — t')dt' = 
{wi ®W2){t)j where ® is the convolution operator. Simi- 
larly, probability of completing N steps at time t is given 

by, 

PNit) ^ {wi® . . .Wi® . . .WN){t). (6) 

Assuming that cancer is an unlikely event (kit <C 1), we 
can use a Taylor expansion in t 

(Tv^TT),!!^^- (7) 

PN{t) is the incidence of cancer; the probability that an 
individual of age t develops cancer in the time interval dt. 



An (invading heterogenous population asymmetric) stem 
cell population will age at a lower rate, k[ = pki. The 
likelihood of incidence of cancer is smaller by a factor of 
. For p— 0.9, and 0.4, and a typical value of = 5 
3J, incidence of cancer is roughly reduced by 40%, and 
99% respectively. 

Fitting the data 

To fit the simulation, we numerically solved Eq|4] us- 
ing the initial condition N{rg/f,t = 0) = A^G'(0.5 ± a), 
where G is a Gaussian distribution with mean 0.5 and 
standard deviation a, and N is a, system-size dependent 
fitting parameter. The other fitting parameter is the pro- 
portionality constant a in EqjSj For Fig.4A, N = 85000 
and a = 1.15. In above scheme, the random nature of 
the lattice (for our choice of randomness -30% of lattice 
constant- mean coordination number is 6 with standard 
deviation 1) biases the average aging rate p so that a ho- 
mogeneous population {(t—0) will have p — 1.015. When 
a hexagonal lattice with no randomness is used, p=l. A 
0.015 offset is added to the theoretical curves to account 
for this. For Fig.7C (top) N = 200000 and a = 1.25. For 
each figure, the four curves shown are fitted simultane- 
ously using the same two fitting parameters. 
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